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Long-range order and low-energy spectrum of diluted 2D quantum AF 
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The problem of diluted two-dimensional (2D) quantum antiferromagnet (AF) on a square lattice 
is studied using spin-wave theory. The influence of impurities on static and dynamic properties is 
investigated and a good agreement with experiments and Monte Carlo (MC) data is found. The 
hydrodynamic description of spin waves breaks down at characteristic wavelengths A > exp( """"^ ), 
X being an impurity concentration, while the order parameter is free from anomalies. We argue that 
this dichotomy originates from strong scattering of the low-energy excitations in 2D. 

PACS numbers; 71.10.-b, 75.10.Jm, 75.10.Nr, 75.40. Gb 



The interest in magnetic properties of high-Tc com- 
pounds has been a major driving force of intensive 
studies of low-dimensional magnetic systems during 
the last decade Q. One of such systems is the 
La2Cui_a;Zn(Mg)a;04, a parental compound doped with 
static vacancies, which shows greater stability of the AF 
order against doping than its mobile hole doped counter- 
part La2_2;Sr2;Cu04 j^J^. It represents a fine example 
of a diluted 2D quantum S = 1/2 AF whose proper- 
ties are the subject of this study. The problem of di- 
luted spin systems has attracted much attention in the 
past and their physics is well understood. The gen- 
eral understanding is that the low-energy excitations of 
these systems are weakly damped spin waves which be- 
long to the infinite cluster |^ and are well defined up 
to the percolation threshold. In this work we show that 
quantum effects and low dimensionality prevent the ex- 
citation spectrum to be defined in these hydrodynamic 
terms at arbitrarily small doping and lead to a para- 
doxical situation where the long-range order is preserved 
but the long- wavelength spectrum is not ballistic. 

The site-diluted AF on a square lattice is described by: 



(1) 



with Pi = 1(0) for the magnetic (non-magnetic) site. The 
pure system at long wavelengths can be described by 
the nonlinear cr-model 0]. Its applicability in the pres- 
ence of a quenched disorder has been questioned Q since 
impurities destroy the Lorentz invariance of the model. 
Some generalizations of the cr-model have nevertheless 
been proposed with parameters modified according to 
MC data |^ and the classical percolation theory Q . 

In this work we study the problem of disorder using 
the i-matrix approach combined with a configurational 
average over the random positions of impurities. It also 
leads to description of the system in terms of an effec- 
tive medium with renormalized parameters. In an ear- 
lier work Harris and Kirkpatrik ^ have shown that in 
2D a non-hydrodynamic term appears in the spin-wave 
self-energy which explicitly violates Lorentz invariance: 
Sk(w) c>c a;/c(ln|a'| — i7r/2). A recent study by Wan et 



al. has confirmed this result. One can show that the 
spectrum is overdamped at wavelengths A ^ exp(7r/4a;) 
which implies the existence of a new lengthscale in the 
system. The absence of a well-defined long-wavelength 
mode also suggests that the true order should be unsta- 
ble under an infinitesimal impurity doping, or that d = 2 
is the upper critical dimension for this problem 0|. 

In this work we study the effect of impurities on the 
spectrum, staggered magnetization, M{x), and Neel tem- 
perature, Tn(x), to clarify the problem of stability of a 
(quasi-)2D AF order. We show that neither M{x) nor 
Tn{x) possesses anomalous contributions which would 
imply an instability. Thus, we have a somewhat paradox- 
ical situation when the spectrum is ill-defined while the 
order parameter is not affected. Such a dichotomy comes 
from the strong infiuence of disorder in the low-energy 
excitations in 2D. The averaging procedure, which effec- 
tively restores translational invariance, fails to recover 
the long-wavelength excitation spectrum of this effective 
medium [ |To| . The low- frequency modes do exist in some 
form but they cannot be classified in terms of an effective 
wave- vector and thus the long- wavelength propagation is 
entirely diffusive. This scheme implies an existence of a 
disorder-induced energy scale which restricts the appli- 
cability of the continuum approach. 

We start with the Hamiltonian in the spin-wave 



approximation H = Ho + "H^^^p + H„„p: 
Ho=^45j5]^k(4"k + /3i/3k) , (2) 

k 

ntp^^ASj J2 e'^'^"'^'^'''[n,kk'4«k' (3) 

;eA,k,k' 

+ <kk' + V2i^, {alpl^, + H.c.)] , 

where the quadratic part of the pure host Hamilto- 
nian TYq is diagonalized using Bogolyubov transforma- 



tion: 
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= Uk^k + Wk/3lk: Sl g/V^S 



UkPl + WkQ:_k, A and B denote the sublattices. 



with 



1, 2wkWk 



-7k/t^k, 7k = (cos(fcj;) -I- 



cos(fcjy))/2, and bare spin- wave frequency ciJk = yl 
All momenta belong to the magnetic Brillouin zone 
kx + ky < IT. Tif = T^rmpi^ ~^ interactions are 
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given by = Vg^ = uiU2+ji-2ViV2+"fiViU2+j2V2'Ui: which is free from the zero-frequency pole. 



V^'^ = Vg"" = V1V2 + 7i_2UiU2 + I1U1V2 + 72U2W1, and 

^A,i,2 = ^s,2,i = "1^2 + 71-2W1W2 + 7i«i^'2 + 72W2M1, 
I runs over the impurity sites. As a first step we solve 
the scattering problem for the single impurity. Since the 
impurity potential is short-ranged one uses the t-matrix 
approach. The advantage of this method is that the sin- 
gle impurity problem can be solved exactly (within the 
spin- wave approximation) |pj| . Therefore, the results of 
this approach become exact in the dilute limit x 0. 

In our case the square lattice site-defect scatters only 
in S-, Px-, Py-, and c?-symmetric channels and one needs 
to solve the t-matrix for each harmonic. The t-matrix 
equations are: 



raoi 
A,i 

pa/3 
^ A,i 



-V. 



A,i 
AA 



A,i ^a'^ A,i 



/3 A,i ' 



A,i '-^a^ A,i ; 



(4) 



where V^'s are the i — s,Px,Py,d components of V's (|^). 
We suppress k and uj in vertex functions rj(k, k', cj), 
bare Green's functions G°(k, cj) = G'^(k, — tj) = {lu — 
Wk + «0}~^, and interactions 14;k'- Products VGT in- 
volve summation over the internal momenta, all ener- 
gies are in units of AS J. An equivalent set of equations 
gives the vertices F^'^^. For the impurity in B sublattice 
TB,i = TA,iia <-> P). Since the impurity potential in each 
partial wave is separable: V^" = Wk'-^k'Uk'^k', 
0fe„Wk'/'fe;Wk'/2, = 7k"k7k'«k' 



sin(fc<^), 7i7 



^ A,s 

ap _ 
A,d 



-tikt^kWk'Wk', v. 



7k«k7k''"k' 



Q/3 

where (t>k^ 



(cos(fca;) — cos(fcj,))/2, Eqs. can be readily solved: 
-yj'j'(k,k')f,(w), where 



r^^(k,k',^ 



fs[Lo) = -l/uj - (1 + uj)p{lo)/[1 - UJ{1 + uj)p{oj)] , 
fp{Lu) - 2/[l +uj + {l- uj){lj^p{uj) - pd{uj))] , (5) 
fd(w) = l/[l + {l-uj)pd{uj)] , 

with p = Ep - u^D and p, = Ep(7p ~ u;p, 

which can be expressed through the complete elliptic in- 
tegrals For the impurity in B sublattice T]f^ — 

— Vgl Ti{—uj). The s-wave scattering (H) reveals a zero- 
frequency mode which originates from the oscillations 
of the fictitious degrees of freedom at the impurity site. 
Roughly speaking, since in the spin-wave approximation 
spins are quantized through bosons even if S' at the im- 
purity site is set to zero there is still ajap left in Sq (for 
discussion see Ref. §). This gives rise to the unphysi- 
cal zero-frequency mode which has to be projected out. 
We do so by introducing magnetic fields at the impu- 
rity sites (similar to Refs. [|l|,|l|l) = H.J^iaU^- 
Within this approach after straightforward algebra 
in the limit 00 one obtains: r^°'g(k, k',a;) — 

-FX-^; (k, k')f. (a>) + AF^: (k, k', 



T.ioj) = -(1 + lo)p{lo)/[1 - ujil + uj)p{uj)] 



(6) 



MkUk' [Wk + ^k' 



AT' 



aP 



-WkUk'l^^k-Wk'+w]. Fs 



-UkWk'[^k + ^k' + H> 

= F^(— w){a ^ /3} as 



AF 

before, p and d waves are not affected by the projection. 

The next step is the averaging procedure which restores 
translational invariance. Assuming random distribution 
of impurities one readily transforms scattering vertices 
into the spin- wave self-energies: Eg"^ (w) = Ei^fk (^)' 
with (co) = x^k-k' [TTl (k, F--; (k, k', 



^^ZH/x = -col[ulf,{u;) + vlTsi-u;)] + 2cOk - , 

Kidiui^) = ^K^'^ H^Pid)i^) + 4f,(,)(-^)] , (7) 
S«^(a;)=xB^f,(.;)+f,(-c.)] , I]^^(c.) = E^£(-^), 

where = -col + (7k K = -(7k = ^k^k^^k, 

^k — ^k^k^'k, B^ = A^uy^Vk- It is interesting that "on 
shell" {lu = (jjk) Eqs. (^ and (^) yield identical Es,k('^k)- 
At low energies oj, cjk self-energies are given by: 

E£"(w) ~ xLOk [p{lo) + 2- 7r/2] ~ xlo + 0{uji,uj^p^) , 

(w) ~ xtOk[piL^) + 7r/2] + Oicokto^p^) , (8) 
with p{uj) ~ (2/7r) ln(|cj|/4) - i , 

which includes contributions from s- and p-wave scatter- 
ing, Ed ~ O(tj^). One can see that besides the "normal" 
softening of the long-wavelength mode and the damp- 
ing proportional to the higher power of k (~ k'^) one 
acquires a non-linear dispersion term with the damp- 
ing 7k /"^k ~ X having only parametric smallness with 
respect to the bare spectrum. Renormalization of the 
real part of the spectrum is dominated by the ln(|ti;|) 
term at low frequencies. A naive "on-shell" pole equation 
'I'k = ijJk + 5iEJJ"(a;k) suggests a vanishing of the spec- 
trum 1^ and an instability of the ground state. Consider- 
ation of the diagrammatic series for the Green's functions 
shown in Fig. |l| with self-energies defined in Eq. (Q) and 
Gi\Lo)^{co-io^-i:^"iu)}-' gi 



ives: 



[GiH^)Gi\-u;)]-^-{K''i^W ' 



Gfiuj) = Ti^^ 



[Gi\co)GlH~co)]-^~i^f{co)r 



(9) 



At fc 3> Wo ^ exp(— 7r/4a;) spectral function reveals 
damped quasiparticle peak at lu ^ k plus some structure 
at low frequencies u loq. At low k ^ luq the quasiparti- 
cle peak merges with this structure and disappears. The 
spectrum is completely overdamped at k uo- Beyond 
this scale no hydrodynamic description is possible. 

A potential concern is whether an anomaly in Eq. (|^) 
can be fixed by the higher order terms. General discus- 
sion of this problem in terms of the susceptibility (Ref. 
H) concluded that the divergency should be unrelated to 
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the expansion in x. One can also show that the summa- 
tion of the Green's function series (^) and dressing of the 
inner hnes in self-energies using self-consistent equation 
on p{uj) and G"^*^ actually enhance the anomaly. 

In a complementary problem of a quasi-2D AF a 
small inter-plane coupling a = J±/J cuts off the log- 
singularity: P3d{^jj) — (I/tt) ln(a/32) -I- at < 
\/2a. Therefore a "safe" range of concentrations x < 
X ^ ln~'(l/a) can be found where the long-wavelength 
quasiparticles are still well defined deep in the 3D region 
of the k-space (A~^ <C \/ct). However, one should be able 
to observe a nonlinearity of the spectrum and an abnor- 
mal damping of the spin waves in the 2D long-wavelength 
region {y/a < A~-^ <C 1) similar to the quasi-lD problem 
For the real materials a ^ 10~^ giving x* ^ 20%. 
Above the concentration x* all the low-energy excitations 
are incoherent because the 2D disorder-induced scale ujq^ 
is shorter than the 3D length l/y/a so the spin waves lose 
coherence before they can propagate in 3D. The same 
consideration applies to the case of small anisotropics in- 
troducing gaps in the spectrum with a modified a — aef / 
accumulating the total effect of the gaps and 3D cou- 
pling. It should be noted that the incoherence comes 
from the averaging procedure which converts the dissi- 
pation of momentum into the dissipation of the energy. 
Therefore, the overdamped excitations should be under- 
stood as diffusive. It is interesting that it requires 2D 
and "strong" disorder to restrict the number of Euclidean 
paths for spin waves and to break down the description 
of the problem in terms of an effective medium. 

Now we proceed with quantities whose expansion in x 
can be shown to be reliable. Staggered magnetization at 
the magnetic site M{x) = Mq - SM (x), Mq ^ S - SX, 

SM{x) = Y^iialau) - 7k(44>]M = (10) 

k 

n{uj)duj p 

where = X^k'^k — 0.1966 is from zero-point fluctua- 
tions, n{uj) = [e"/"^ — 1]"^ is the Bose occupation number, 
index R denotes retarded. At T = n{uo) = —9{—lo) 
and integrals in (|l^) with G"^*^ from (^) can be taken 
numerically. It can be shown [p^ that the expansion 
in X for G'^'^ can be performed before the integration 
GO -h G^I]""G0, G"^ ~ G°S""G0 and that 
all integrals for the linear in x term are convergent 
giving 5M{x) = 0.209(8) • x + 0{x'^). For 5=1/2 
M{x)/Mo ~ I - 0.691(5) • X. This linear slope together 
with the results of calculations using G""' from Eq. (|), 
MC ^ and NMR data (l7| are shown in Fig. |. Note 
that Fig. I and Eq. (|l^) are showing the reduction of the 
magnetic moment by the quantum fluctuations induced 
by impurities. In order to extract the same quantity from 
the MC data ||l^, which are averaged over all sites, one 
needs to deduct a probability of finding a spin-occupied 



site from them. One can see a very good agreement of 
our results with numerical data up to high concentra- 
tions. The oxidation of the crystals can be the reason of 
a faster decrease of M{x) in NMR data. The absolute 
value of SM{x) is independent of S in the linear spin- 
wave approximation similar to the quantum reduction of 
S by ^A. We plot SM{x) in Fig.|^ (right axes) in order to 
emphasize the agreement with the MC data for 5* = 1/2 
and S — 1, which show only weak S'-dependence. 

The quantity which one would expect to exhibit a dif- 
ferent behavior is T^r for quasi-2D problem. The ther- 
mal corrections to the staggered magnetization SM'^ ^ 
rin(T/w)L^^ yield a finite Tn ~ ln(l/a)-i for a pure 
system. Since the original spectrum is bent from the 
linear form this should manifest itself in the finite tem- 
perature part of Eq. (p^ ) as a stronger divergency. One 
indeed finds such contributions in Eq.(^o|) (~ a;Tln^ \uj\) 
already in a perturbative limit (j;ln|w| ^ 1). However, 
these anomalous terms from diagonal and off-diagonal 
parts cancel each other. From the mean-field equation 
M'^'^ (x) — after some algebra one obtains (at a 0): 
TN{x)/TNiO) ~ l-x^ with A = 7r-2/7r-|-0.209(8)/Afo. 
For S = 1/2, A = 3.196(5) |l|l, which fits very weU 
experimental data ]T9| , p0t , see Fig.|[ 

Thus the order is preserved up to a high x value [p^|J^] 
while the long-wavelength spectrum is not well defined 
at any x. Experimentally one would need to probe di- 
rectly the spin-wave spectrum close to the AF ordering 
vector and look for anomalous broadening and nonlin- 
earity. The actual overdamped part of the spectrum is 
associated with a small energy scale (and small momenta, 
e.g. a.t X = 0.2, feg ~ O.OItt/o) and it is not clear if it is 
observable within the range of doping where conclusions 
of this work are reliable. The observed deviation from 
the simple exponential behavior of the correlation length 
^(T, a;) v.s. 1/T ^ might be related to the character of 
the low-energy spectrum we discuss. 

Another aspect relates our problem to the other prob- 
lems of disorder in 2D. We observe that some weight from 
every mode is transferred to the low energies w ~ 
which can be understood as some sort of localization, 
though the localization criteria for this problem is un- 
clear. Altogether these low-energy states result in a peak 
in the density of states. From x-expanded form of the 
Green's function G"" ~ G° -f Gjii;°"G° one readily gets 
N{uj) = 2uj/tt + xC + O{xoj\n\uj\), which implies a fi- 
nite density of states at w = 0. With G"" from Eq. (||) 
one obtains N{uj) « 2Lu/7r + xC/[{l + 4:x\n\Lj\/7r)^ + 'ix^] 
which is zero strictly at w = but has a sharp peak ~ 1/x 
a,t oj = ojQ. This is very close to the dispute over N{uj) for 
the certain types of disorder in 2D systems with linear ex- 
citation spectrum [|l^. Note that in our case there is no 
"weak" limit for disorder since the perturbation (^ JSz) 
is of the same strength as the pure coupling [|l| . 

In conclusion, we have provided evidence for the sta- 
bility of the long-range order in the 2D and quasi-2D AF 
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doped with static vacancies. At the same time the long- 
wavelength excitation spectrum is shown to be diffusive 
for any value of doping, restricting a hydrodynamic ap- 
proach to the problem from low energies. 
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FIG. 1. Belyaev diagram series for the diagonal, G"" , and 
off-diagonal, G"'', Green's functions. Self-energies (cir- 
cle) and E"'' (square) are the configurational averages of T"" 
and r"'', respectively. 
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FIG. 2. Left axes: reduction of the local magnetic mo- 
ment by the quantum fluctuations induced by impurities 
M{x)/M{0) from NMR (o, Ref. [0) and MC (•, Ref. @) 

studies. Lines show our results from Eq^ (p^. Right axes: 

the absolute value of SM{x) from Eq. ( |lo| ) (lines) and MC 
data, Ref. (li , for 5* = 1/2 (o) and S = 1 (<). 
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FIG. 3. Neel temperature v.s. x from fiSH and susceptibil- 
ity measurements Refs. p9|j20) for the La2Cui_iZn(Mg)a;04, 
and the linear slope from our theory. 
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